A unified origin for the 3D magnetism and superconductivity in Na^CoC^ 
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We analyze the origin of the three dimensional (3D) magnetism observed in nonhydrated Na-rich 
Na^Co02 within an itinerant spin picture using a 3D Hubbard model. The origin is identified as the 
3D nesting between the inner and outer portions of the Fermi surface, which arise due to the local 
minimum structure of the a\ CJ band at the F-A line. The calculated spin wave dispersion strikingly 
resembles the neutron scattering result. We argue that this 3D magnetism and the spin fluctuations 
responsible for superconductivity in the hydrated systems share essentially the same origin. 



The discovery of superconductivitWSC) in NazCoCVy 
H2O has attracted much attention. [l| Although the ex- 
periments are still somewhat controversial, we briefly 
summarize the present understanding, (i) Up to now, 
all the angle resolved photoemission (ARPES) measure- 
ments E3 [1 show the absence of the e' g hole pock- 
ets predicted in the first principles calculation. 6] Thus 
the remaining a\ g seems to be the relevant band, (ii) 
In the bilayer hydrated (BLH) SC samples, there is an 
enhancement in QVT) -1 at low temperatures at the Co 
site 0, S H [H K3 and at the O site, d E3 (Tf 1 is 
the spin-lattice relaxation rate) indicating the presence 
of spin fluctuations (SF) located away from the Bril- 
louin zone (BZ) edge.|9l Il2| while such an enhancement 
is not seen in non-SC monolayer hydrated and nonhy- 
drated Na-poor samples. On the other hand, the Knight 
shift stays nearly constant at low temperatures in the 
SC samples, 0, 0; 113 which means that the SF is not 
purely ferromagnetic. The bottom line is that SF that 
is neither purely ferromagnetic nor antiferromagnetic is 
strongly related to SC. (iii) Unconventional SC gap is 
suggested from the absence of the coherencepeak in T-T 1 
as well as the power law decay below T c .0, H, EL pji 
Several experiments suggest spin singlet pairing. |l4 llq 
and moreover, the effect of the impurities on T c is found 
to be small, (jjjllli suggesting a s-wave like gap. 

In a previous paper, as a solution for understanding 
these experimental results, we have proposed a mech- 
anism for an unconventional s-wave pairing, in which 
the local minimum structure (LMS) of the a± g band at 
the r point plays an important role.|l9j For appropri- 
ate band fillings, LMS results in disconnected inner and 
outer Fermi surfaces (FS), whose partial nesting induces 
SF at wave vectors that bridge the two FS (see FigH3 
left panel). This incommensurate SF can give rise to an 
"extended s-wave" pairing in which the gap changes sign 
between the two FS, but not within each FS. A recent 
multiorbital analysis also shows that such a pairin g ca n 
occur when the inner and outer a\ g FS are present. |20j 

The purpose of this paper is to investigate the relation 
between this SC and the 3D spin-density-wave(SDW)- 



like, metallic magnetism observed in the nonhydrated 
Na-rich systems, |2l|, |22j, |23j, \2A\ whose magnetic struc- 
ture has been revealed by neutron scattering experiments 
to have in-plane ferro mag netic and out-of-plane antifer- 
romagnetic character. |2j,|2fJ Neutron scattering exper- 
iments have further obtained the spin wave dispersion 
of this 3D magnetism. Analysis of this dispersion based 
on the localized spin picture have found that the out-of- 
plane antiferromagnetic coupling is of the order of the 
in-plane ferromagnetic coupling, 0, EE E3 despite the 
strong 2D nature of the material. Moreover, Curie- Weiss 
temperature is expected to be positive from the evalu- 
ated co upling constants . |25j while it is actually negative. 
[3 Ulllll El EHH HH A theory based on a spin-orbital 
polaron picture has been proposed for this puzzle. |3(ij | 

In the present study, since the system remains metal- 
lic in the magnetically ordered state, we take an itiner- 
ant spin viewpoint. From the consistency between the 
experiments and the calculated results, we propose that 
the origin of the 3D magnetism is essentially the same 
with the SF responsible for SC : the nesting between in- 
ner and outer portions of the Fermi surface that arises 
due to the local minimum of the a\ g band. [3jJ 

There exist two C0O2 layers within a unit cell due 
to the alternation of the oxygen arrangement, but nei- 
ther the effective hopping integrals nor the on-site en- 
ergy alter along the c axis, so the BZ folding in the c 
direction occurs without inducing a gap in the band. 
Thus, in a 3D effective theory in which O and Na site 
degrees of freedom are integrated out, the unit cell con- 
tains only one layer, which results in an unfolded BZ 
shown in the upper right of FigQ Taking into account 
only the a\ g band, we consider a single band Hubbard 

model H = E tJ E CT ^ c L c j<x + C/Ei n *T n a> 011 a 3D 
triangular lattice, where the 1st, 2nd and 3rd neigh- 
bor in-plane hopping integrals t±, t%, £3 and the out-of- 
plane nearest neighbor hopping t z (Fig^ upper left) are 
chosen so as to roughly reproduce the a\ g portions of 
the band obtained in the first principles calculation. 
Throughout the paper, we take t 2 = —0.35, ts = —0.07, 
and t z = —0.15 (unless otherwise noted) in units of 



2 






'.LMS 



km r 



FIG. 1: (color online) Upper left: the 3D model in the present 
study. Upper right: the Brillouin zone. The wave numbers in 
the parenthesis are those in the unfolded BZ scheme. Lower 
panel: the band structure of the present model with U — 0. 
k z = 0, 7r/2,7r are k z in the unfolded BZ scheme. 



t\ = 1, which corresponds to about (or slightly less 
than) O.leV according to ARPES measurements. |4| The 
band (with U = 0) for this choice of parameter values is 
shown in Fig^ The band filling is n =number of elec- 
trons/site, and it is related to the actual Na content x 
by n — 1 = x (provided that the e' g bands are fully occu- 
pied). We calculate the spin susceptibility using the fluc- 
tuation exchange (FLEX) approximation^ as Xs(q) = 
Xirr(<7)/[1 - Uxin{q)], where q = (q,iu n ). Here, the irre- 
ducible susceptibility is Xin (<z) = — ■h Y^,k + l)G(k) 
(iV:number of fc-point meshes), where G is the renormal- 
ized Green's function self-consistently obtained from the 
Dyson's equation, in which the self energy is calculated 
using G and Xirr- We take up to N = 64 x 64 x 64 fc-point 
meshes and up to 16384 Matsubara frequencies. 

In order to analyze the magnetically ordered state with 
a wave vector Q at which the FLEX spin susceptibility 
is maximized, we consider a mean field Hamiltonian 
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FIG. 2: (color online) The vertical cross section (k y = 0) of 
the FS for various n with [7 = 0. Here we use the unfolded 
BZ scheme, but for n — 1.60 we also show the FS in the usual 
folded BZ scheme by dashed lines. The unfolded BZ edge in 
the k z direction is denoted as "^yj" to stress that the unit 
cell is halved in the c direction. The finite thickness is due to 
a finite range of energy [Ef ± 0.03£i) taken in obtaining the 
FS. The thickness thus represents the density of states or the 
effective mass. The arrows present the diameter of the FS. 
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FIG. 3: (color online) Similar as in Fig|3 except t2 = t$ = 0. 
The band is shown in the inset. 



where s(k) is the bare band dispersion, cj e+m Q a = 
2 u a,o-,™(fc)7fc a a is a unitary transformation that 
diagonalizes the Hamiltonian to obtain the two bands 
E\ a and E^a, and (tiq a ) is self-consistently determined 
in the usual procedure. We then calculate the irreducible 
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FIG. 4: (color online) Upper left: The band filling dependence 
of Tsdw obtained by FLEX. The inset shows the T depen- 
dence of Xs(0,0) -1 for n = 1.82 and n — 1.7. Upper right: 
The FLEX x= [q, 0) for n = 1.82, U = 6, and T = 0.05. Lower 
left: the nesting vector of the FS for n — 1.82 is shown. Lower 
right: The FS in the ordered state (T = 0.015) obtained by 
the mean field calculation is shown. Although the unit cell 
is "truely" doubled in this case (tt/2 < \k z \ < n part should 
actually be drawn repeatedly in < \k z \ < n/2 region), we 
continue to use the unfolded scheme for a clear comparison 
with the FS in the non-ordered state in the left. 



susceptibility matrix by 

P rn,n,a,a' 

^ u a,a : m-\-n-\-j (jP H~ &)^a' ,0' .m-{-n {jPj^ot' ,a' ,rn (jp) 

x /(iW(p))-/(£ a , CT (p + fc)) 
E a , a {k + p) -E a i^iip) +ui + ie' 

where f(E) is the Fermi distribution function. The 
spin wave dispersion can be determined by the condi- 
tion that the real part of the eigenvalue of the ma- 
trix 1 — Uxl (q,u)) equals zero. We take up to N — 
192 x 192 x 96 fc-point meshes here. 

We first show in FigElthe evolution of the 3D FS (with 
U = 0) as the band filling is increased. Here the vertical 
cross section of the FS is given in the unfolded BZ. The 
FS is a 2D cylinder for small band filling, that is, for low 
Na content. As the band filling is increased, an inner FS 
appears around the r point. As in the purely 2D(~ BLH) 
case (see FigEJor ref ll9|) . this FS is disconnected from the 
outer 2D cylindrical FS, but the inner FS is 3D in the 
NH case. 33] For larger band fillings, the two FS become 
connected (around n = 1.73), and then finally for higher 
band fillings it becomes a single 3D FS. The appearance 
of the 3D FS despite the small t z is a consequence of LMS 
of the band along the T-A line. In fact, if we take a band 



that does not have LMS, i.e., for = £3 = but with the 
same t z (— —0.15), the FS remains to be 2D (unless the 
band filling becomes very close to 2) as shown in FigOU 
The 3D evolution of the FS is consistent with the ARPES 
observations, where the in-plane kp barely decreases with 
the increase of the Na content for x > 0.6. 0,13 In fact, 
the decrease of the FS diameter in Fig [2 is very small 
compared to that in FigOU 

In the inset of the upper left of Fig0] , FLEX re- 
sults of X7 1 (Q = 0,zw„ = 0) are plotted as functions 
of T for U = 6 and n = 1.82 or n = 1.7. The Curie- 
Weiss behavior at high temperatures with a negative 
Curie- Weiss temperature is consistent with the exper- 
iments, d m in m m m m m in the upper 

left of FigOJ Tsdw, defined as the temperature where 
the maximum value of U\irr reaches ~ 1, is plotted 
as a function of the band filling. Tsdw exists only 
in the regime 1.75 < n < 1.9, which is consistent 
with the fact that the magnetic ordering is not ob- 
served experimentally for x < O.75.0,|2^| The maximum 
?sdw corres pon ds to about 20K, in agreement with the 
experiments. |2ll E2I E3 | Around n ~ 1.8, the ordering 
occurs with Q = (0, 0, it), as shown in the upper right 
of FigQ] for n — 1.82, which corresponds to in-plane 
ferromagnetic, out-of-plane antiferromagnetic spin struc- 
ture, as found experimentally. |2fl Ef^ The origin of this 
spin correlation can be found in the shape of the 3D FS. 
Namely, the FS at n ~ 1.8 is partially nested with a 
nesting vector close to (0,0, n), as can be seen in FigQ] 
(lower left panel). In fact, this nesting occurs between 
portions of the FS that can be considered as remnants of 
the inner and outer FS, which are disconnected for lower 
band fillings. |3fl| On the other hand, as the band filling 
comes close to n — 1.9 {x = 0.9), the FS becomes too 
small for (0, 0, tt) nesting, and the ordering wave number 
becomes incommensurate in the c direction, which is also 
consistent with the finding in refl22l 
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FIG. 5: (color online) The calculated in-plane (left) and out- 
of-plane (right) spin wave dispersion for n = 1.82, [7 = 4, and 
T = 0.015. Here the wave vectors are presented in units of 
the reciprocal lattice primitive vector in the usual folded BZ 
for clear comparison with Fig.4 of refl25l. 
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ened if the two dimensionality is increased by hydration, 
resulting in enhanced SF. As we have mentioned in the 
beginning, our proposal is that this SF is responsible for 
SC.|l9|| Thus, in our view, the 3D magnetism and the 
incommensurate SF that gives rise to SC share the same 
origin as shown in Fig|S] the nesting between the inner 
and the outer (connected or disconnected) FS that orig- 
inate from the local minimum of the a\ g band. 

Numerical calculation were performed at the facilities 
of the Supercomputer Center, ISSP, University of Tokyo. 
This study has been supported by Grants-in-Aid for Sci- 
entific Research from the Ministry of Education, Culture, 
Sports, Science and Technology of Japan, and from the 
Japan Society for the Promotion of Science. 



FIG. 6: (color online) The FS nesting relevant for SC and the 
3D magnetism. 



We now move on to the mean field results in the 
Q = (0, 0, tt) SDW ordered state. We focus on n = 1.82 
and take U smaller than U = 6 adopted in FLEX since 
the mean field approach tends to overestimate the ten- 
dency toward ordering. In fact, taking U = 4 gives a 
self-consistently determined (tiq^) value of 0.075 at low 
T, corresponding to a magnetic moment of 0.15/xb per 
site (Co atom) . in rough agreement with a /iSR estima- 
tion 0.18/xs.|22| In Fig0|(lower right panel), we show the 
calculation result of the FS in the ordered state. The por- 
tions of the FS around k z — ±tt/2 disappear, but other 
portions remain to result in a FS with a strong 2D charac- 
ter. This is consistent with the fact that system remains 
metallic in the magnetically ordered state, [jl) Note also 
that the disappearing portions of the FS are thick, i.e., 
they have heavy mass, which is consistent with what has 
been suggested in reffSll from the enhancement of the 
mobility in the magnetically ordered state. 

In Fig |51 we show the calculated spin wave dispersion. 
The overall feature as well as the in-plane to out-of-plane 
ratio w(G.l, 0.1, 0)/w(0, 0, 0.5) ~ 2.1 strikingly resembles 
the neutron scattering result of refl25l for x — 0.82. The 
energy scale (w(0, 0,0.5) ~ O.lti corresponds to ~ 10 
meV) is also close to the experimental result. Here again, 
our results show that the 3D magnetism is well under- 
stood within the present itinerant spin picture. |36| 

Finally, let us discuss the relation between this 3D 
magnetism and SC. As we have seen, the SF in the non- 
hydrated systems becomes weak for small band fillings. 
This is because the nesting between the 2D outer and 
the 3D inner FS is not good. However, the nesting be- 
tween the inner and the outer FS[37j can be restrength- 
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